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The  interaction  of  a  vortex  pair  with  a  premixed  flame  serves  as  an  important  prototype  for  premixed 
turbulent  combustion.  In  this  study,  we  investigate  the  interaction  of  a  counter-rotating  vortex  pair  with 
an  initially  flat  premixed  methane  flame.  We  compare  computational  solutions  obtained  using  a  time- 
dependent,  two-dimensional  adaptive  low  Mach  number  combustion  algorithm  that  incorporates  GRI- 
Mech  1.2  for  the  chemistry,  thermodynamics,  and  transport  of  the  chemical  species.  We  discuss  the  me¬ 
chanical  nature  of  the  flame-vortex  interaction  and  the  features  of  the  interaction  that  are  strongly  affected 
by  fuel  equivalence  ratio,  f.  We  find  that  the  circulation  around  the  vortex  scours  gas  from  the  preheat 
zone  in  front  of  the  flame,  making  the  interaction  sensitive  to  equivalence  ratio.  For  cases  with  r/j  *=  1 ,  the 
peak  mole  fraction  of  CH  across  the  flame  is  relatively  insensitive  to  the  vortex,  whereas  for  richer  flames 
we  observe  a  substantial  and  rapid  decline  in  the  peak  CH  mole  fraction,  commencing  early  in  the  flame- 
vortex  interaction.  The  peak  concentration  of  HCO  is  found  to  correlate,  in  both  space  and  time,  with  the 
peak  heat  release  across  a  broad  range  of  equivalence  ratios.  The  model  also  predicts  a  measurable  increase 
in  C2H2  as  a  result  of  interaction  with  the  vortex  and  a  marked  increase  in  the  low-temperature  chemistry 
activity. 


Introduction 

Turbulence  affects  the  process  of  methane  com¬ 
bustion  through  a  wide  variety  of  mechanisms.  Tra¬ 
ditional  approaches,  based  on  asymptotic  analysis, 
show  that  velocity-induced  tangential  strain  at  the 
flame  surface  can  dramatically  enhance  or  suppress 
combustion  activity  in  the  flame  zone  depending  on 
Lewis  number.  These  effects  have  been  studied  us¬ 
ing  results  from  steady  flat-flame  counterflow  exper¬ 
iments  (see  Refs.  [1-5]).  As  discussed  in  the  review 
by  Peters  [6],  this  type  of  information  can  be  readily 
incorporated  into  engineering  models  through  the 
flamelet  concept. 

Research  aimed  at  understanding  transient  effects 
in  turbulent  flows  has  led  investigators  to  study  the 
interaction  of  a  single  vortical  structure  with  a  pre¬ 
mixed  flame.  These  flows  offer  advantages  of  being 
more  repeatable  and  more  amenable  to  analysis  and 
experimental  diagnostics  than  fully  developed  tur¬ 
bulent  flow.  Such  studies  typically  consider  either  a 
planar  vortex  pair  (e.g.,  see  Refs.  [7-11])  or  an  axi- 
symmetric  toroidal  vortex  (e.g.,  see  Refs.  [12-19]). 
Vortex-driven  flame  quenching  has  been  studied 
computationally  using  both  single-step  [20]  and  two- 
step  [21]  reaction  mechanisms.  The  performance  of 


certain  flame  markers  in  the  presence  of  time-de- 
pendent  vortical  perturbations  has  been  assessed  in 
simulations  using  skeletal  Ci  methane  chemistry 
[7,22]  and  more  complete  description  of  Ci  and  C2 
chemistry  [23].  Hilka  et  al.  [24]  assess  the  ability  of 
a  4-step  mechanism  to  predict  combustion  heat  re¬ 
lease  in  this  configuration,  using  a  more  detailed  52- 
step  mechanism  as  a  baseline. 

While  most  of  the  experimental  data  on  methane 
flame  quenching  and  combustion  markers  focus  on 
lean  flames,  Nguyen  and  Paul  [8]  considered  a  rich 
methane  flame.  Their  results  are  particularly  inter¬ 
esting  for  several  reasons.  First,  they  predict  dra¬ 
matic  changes  in  species  concentrations,  notably  CH 
and  OH,  that  are  not  predicted  by  counterflow  flame 
simulations  [2],  This  suggests  that  a  flamelet  model 
based  solely  on  flat-flame  data  may  be  inadequate 
for  transient  flow  field  interactions.  Second,  the  lim¬ 
ited  number  of  computational  studies  available  for 
modeling  rich,  multidimensional  vortex-flame  inter¬ 
actions  yield  results  that  are  somewhat  inconsistent 
with  the  experiment.  For  example,  Najm  et  al.  [23] 
computed  a  flow  problem  similar  to  the  Nguyen  and 
Paul  experiment.  They  used  the  GRI-Mech  1.2  [25] 
reaction  mechanism  and  a  smaller,  faster  vortex  with 
less  nitrogen  dilution;  however,  their  simulations 
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LAMINAR  PREMIXED  FLAMES — Vortex  Flame  Interactions 


failed  to  reproduce  several  key  features  observed  in 
the  experiment. 

In  this  paper,  we  study  in  detail  the  interaction  of 
a  vortex  with  premixed  N2-diluted  methane  flames 
that  have  equivalence  ratios  in  the  range  <j)  =  0.76- 
1.28.  The  vortex  we  impose  matches  the  strength, 
core  separation,  and  speed  of  that  observed  in  the 
Nguyen  and  Paul  experiment.  Although  our  config¬ 
uration  is  close  to  the  actual  experimental  condi¬ 
tions,  our  initial  premixed  flame  is  flat  and  oriented 
normal  to  the  inlet  flow.  The  simulations  therefore 
lack  a  small  lateral  velocity  component  associated 
with  the  V-flame  configuration  of  the  related  exper¬ 
iment.  We  cany  out  these  studies  by  numerical 
simulation  using  a  parallel,  adaptive  low  Mach  num¬ 
ber  combustion  algorithm  developed  by  Day  and 
Bell  [26],  All  of  the  computations  use  chemistry, 
transport,  and  thermodynamic  databases  specified  in 
GRI-Mech  1.2.  Our  results  match  the  CH  behavior 
reported  in  Nguyen  and  Paul’s  experiment.  In  ad¬ 
dition,  we  observe  that  the  flames  response  to  the 
vortex  is  a  strong  function  of  the  inlet  flow  equiva¬ 
lence  ratio.  We  find  that  mole  fraction  of  HCO  cor¬ 
relates,  in  both  position  and  time,  with  the  heat  re¬ 
lease  of  the  methane  flame.  Finally,  we  find  that  in 
spite  of  the  agreement  in  CH  profiles,  our  simula¬ 
tions  do  not  reproduce  the  experimentally  observed 
behavior  of  OH. 


Numerical  Model 

Our  computational  approach  uses  a  hierarchical 
adaptive  mesh  refinement  (AMR)  algorithm  based 
on  an  approximate  projection  formulation  for  inte¬ 
grating  the  momentum  equations  [27]  and  subse¬ 
quently  extended  to  low  Mach  number  flows  [28]. 
The  projection  algorithm  is  coupled  to  conservation 
equations  for  chemical  species  and  enthalpy.  The 
single-grid  algorithm  is  implemented  in  a  structured 
uniform  grid  setting  and  incorporated  into  an  AMR 
framework  that  employs  a  recursive  time-stepping 
procedure  over  refinement  levels.  We  sketch  the 
model  and  numerical  implementation  below;  the 
reader  is  referred  to  Ref  .[26]  for  details. 

Our  model  is  based  on  a  conservative  form  of  the 
low  Mach  number  combustion  model  introduced  by 
Rehm  and  Baum  [29],  subsequently  derived  from 
low  Mach  number  asymptotic  analysis  by  Majda  and 
Sethian  [30],  We  consider  a  gaseous  mixture,  ignor¬ 
ing  Soret  and  Dufour  effects,  body  forces,  and  ra¬ 
diative  heat  transfer,  and  assume  a  mixture  model 
for  species  diffusion  [31,32],  For  an  unconfined  do¬ 
main,  we  have 

+  V-pUU  =  -Vn  +  V-T  (1) 


SpYm  + 
at 


V  ■  UpY, 


V-pVmVYm  -  com  (2) 


^  +  V  •  Uph  =  V  •  —  Vh 

Ul/  c-'p,mix 

+  2  V-hJpDm  -  — )  VYm  (3) 
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where  p  is  the  density,  U  is  the  velocity,  7r  is  the 
pressure  variation  from  a  uniform  ambient  pressure, 
Ym  is  the  mass  fraction  of  species  in,  h  is  the  mass- 
weighted  enthalpy  of  the  gas  mixture,  T  is  the  tem¬ 
perature,  and  wm  is  the  net  destruction  rate  for  pYm 
due  to  chemical  reactions.  Also,  X  is  the  thermal  con¬ 
ductivity,  t  is  the  stress  tensor,  cp  mix  is  the  specific 
heat  of  the  mixture,  and  hm(T)  and  T)m  are  the  en¬ 
thalpy  and  species  mixture-averaged  diffusion  coef¬ 
ficients  of  species  m,  respectively.  These  evolution 
equations  are  supplemented  by  an  equation  of  state 
for  a  perfect  gas  mixture: 

Po  =  pRm*T  =  p'RT'Z^-  (4) 

m  */Vrn 

where  Wm  is  the  molecular  weight  of  species  m. 

In  the  low  Mach  number  limit,  the  thermody¬ 
namic  pressure  given  by  equation  4  remains  approx¬ 
imately  constant  as  the  flow  evolves.  Differentiating 
the  equation  of  state  in  the  frame  of  the  fluid,  and 
using  the  conservation  equations  to  replace  advec- 
tive  derivatives,  we  obtain  an  elliptic  constraint  on 
the  evolving  velocity  field: 

v-h  =  +  2  />Dmvrm-vftJ 

pCpT  \  m  / 

+  -  2  Try  V  ■  pDmVYm 
P  m  Wm 


where 

W  =  (SmYJWm)-1 

and 


r  —  T  Y  rib  /rlT 

°p,mix  mLU  lm'  U'1 

The  single-grid  scheme  that  forms  the  basis  for 
our  adaptive  algorithm  combines  a  symmetric  op¬ 
erator-split  coupling  of  chemistry  and  diffusion  pro¬ 
cesses  with  a  projection  method  for  incorporating 
the  velocity  divergence  constraint.  First,  equations 
1-3  are  advanced  in  time  using  a  second-order  Go¬ 
dunov  scheme  for  convective  terms  and  a  time-cen¬ 
tered  Crank-Nicolson  discretization  for  diffusion. 
Because  the  transport  coefficients  depend  on  both 
temperature  and  composition,  we  adopt  a  sequen¬ 
tial,  predictor-corrector  scheme  to  guarantee  sec¬ 
ond-order  treatment  of  nonlinear  diffusion  effects. 
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Fig.  1.  Schematic  of  the  premixed  methane  flame-vortex 
problem.  The  shaded  line  represents  the  position  of  the 
flame,  and  the  swirl  lines  represent  the  initial  vortex.  We 
impose  symmetry  along  the  sides  of  the  0.8  X  4.0  cm  do¬ 
main  to  avoid  modeling  the  dotted  region. 

The  chemistry  is  advanced  using  time-implicit  back¬ 
ward  differentiation  methods  [33],  The  implicit  dif¬ 
fusion  and  chemistry  components  of  the  algorithm 
are  time  split  symmetrically  to  ensure  that  the  com¬ 
posite  algorithm  remains  second  order.  The  velocity 
field  resulting  from  the  advection/diffusion/chem- 
istry  step  is  then  decomposed  using  a  density- 
weighted  approximate  projection.  Numerically,  the 
projection  step  requires  the  solution  of  a  variable- 
coefficient  linear  elliptic  equation.  The  component 
satisfying  the  constraint,  equation  5,  updates  the  ve¬ 
locity  field,  and  the  remainder  updates  the  pertur- 
bational  pressure.  The  explicit  procedure  for  the 
treatment  of  advection  terms  necessitates  a  Courant- 
Friedrichs-Lewy-type  (CFL-type)  time-step  restric¬ 
tion.  Since  the  advective  time  scale  is  typically  larger 
than  the  fastest  time  scales  associated  with  the 
chemical  kinetics,  this  does  not  appear  to  be  a  seri¬ 
ous  disadvantage  for  time-dependent  simulations. 

We  note  that  because  of  the  fractional-step  nature 
of  the  overall  approach,  it  is  impossible  to  satisfy  the 
equation  of  state  while  conserving  species  and  en¬ 
thalpy  discretely.  We  add  a  correction  to  equation  5: 

/  SP_ 

TPamb  A 

where  dp  =  p0  —  pamb>  Pamb 's  d  ie  ambien  t  pressure. 
At  is  the  numerical  time  step,  and  y  is  the  ratio  of 
specific  heats.  This  forcing  term  prevents  the  algo¬ 
rithm  from  drifting  an  appreciable  amount  off  the 
equation  of  state.  The  details  of  this  correction  and 
its  effects  are  discussed  in  Ref.  [26];  however,  we 
note  that  in  the  calculations  presented  here  the  ther¬ 
modynamic  pressure  deviates  no  more  than  0.5% 
from  the  ambient  atmospheric  pressure. 


The  extension  of  the  above  algorithm  to  adaptive 
mesh  refinement  is  based  on  a  hierarchical  refine¬ 
ment  strategy.  Our  methodology  uses  a  system  of 
overlaid  grids  with  successively  finer  spacing  in  time 
and  space.  Fine  grids  are  formed  by  uniformly  di¬ 
viding  coarse  cells  in  each  direction.  Increasingly 
finer  levels,  each  consisting  of  a  union  of  rectangular 
grid  patches,  overlay  coarser  grid  levels  until  the  so¬ 
lution  is  adequately  resolved.  An  error-estimation 
procedure  identifies  where  refinement  is  needed, 
and  grid-generation  procedures  dynamically  create 
or  remove  rectangular  fine  grid  patches  as  require¬ 
ments  change. 

The  complete  adaptive  algorithm  has  a  number  of 
desirable  properties.  The  overall  method  is  second- 
order  accurate  in  space  and  time  and  discretely  con¬ 
serves  species  mass  and  enthalpy.  Furthermore,  the 
algorithm  satisfies  a  free-stream  preservation  prop¬ 
erty  such  that  non-reacting  isothermal  flow  remains 
numerically  isothermal  during  species  transport. 
The  methodology  has  been  implemented  for  distrib¬ 
uted  memory  parallel  processors  using  the  BoxLib 
class  libraries  described  by  Rendleman  et  al.  [34], 

Results 

We  conducted  a  detailed  study  of  a  counter-rotat¬ 
ing  pair  of  vortices  interacting  with  N2-diluted  pre¬ 
mixed  methane-air  flames.  We  focus  the  discussion 
on  two  aspects  of  the  results.  First,  we  are  interested 
in  characterizing  the  mechanical  nature  of  the  flame- 
vortex  interaction.  Second,  we  want  to  understand 
how  this  behavior  affects  the  combustion  process  as 
a  function  of  flame  stoichiometry. 

The  computational  study  approximates  the  exper¬ 
imental  conditions  of  Nguyen  and  Paul.  Following 
Najm  et  al.,  we  approximate  the  V-flame  experiment 
as  a  flat  premixed  flame  oriented  normal  to  the  inlet 
flow  and  superimpose  a  velocity  field  due  to  a  peri¬ 
odic  array  of  counter-rotating  vortex  pairs  with  Gaus¬ 
sian  cores  <5V  =  2.25  mm  wide  and  centers  rv  =  0.25 
cm  apart.  The  vortex  pair  propagates  upward  with  a 
self-induced  velocity  of  110  cm/s.  These  parameters 
produce  a  vortex  pair  with  the  same  width  and  prop¬ 
agation  speed  as  reported  for  the  vortex  pair  in  the 
Nguyen  and  Paul  experiment.  We  numerically 
evolve  only  the  left  vortex,  as  indicated  in  Fig.  1 .  Our 
base  computational  domain  is  i',|  =  8  mm  wide  and 
ijd  —  4  cm  high,  with  the  boundary  along  the  right 
side  corresponding  to  the  vortex  centerline.  The  pe¬ 
riodic  configuration  results  from  symmetry  bound¬ 
ary  conditions  imposed  on  the  sides  of  the  domain. 
Reactants  flow  in  at  the  bottom,  and  combustion 
products  exit  through  the  top. 

We  position  the  initial  premixed  flame  so  that  the 
peak  mole  fraction  of  HCO  occurs  at  iff  —  2  cm  (see 
Fig.  1).  Cell-centered  values  are  initialized  using  cell 
averages  of  a  refined  PREMIX  [35]  solution.  A  vor¬ 
tex  is  superimposed  on  this  flow  field,  centered  half¬ 
way  between  the  inlet  and  flame.  We  cover  the  en¬ 
tire  domain  with  a  uniform  coarse  mesh  (Ax  =  250 
/mi)  and  then  place  up  to  two  refined  grid  levels, 
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Fig.  2.  CH  mole  fraction.  In  each 
column,  the  three  frames  corre¬ 
spond  approximately  to  4.3,  6.3,  and 
7.9  ms  simulation  time  from  the  ini¬ 
tial  conditions  (see  Fig.  1).  The  color 
scale  in  each  column  is  normalized 
to  peak  values  ofXCH  =  (0.23,  7.76, 
10.8,  3.66,  0.30)  X  Hr7,  respec¬ 
tively.  The  frames  represent  8  X  9.5 
mm  subregions  of  the  domain  near 
the  flame  surface. 


each  with  Ax  decreased  by  a  factor  of  2,  to  resolve 
the  flame  zone.  The  resulting  resolution  is  equiva¬ 
lent  to  that  of  a  128  X  640  uniform  grid  calculation 
with  Ax  —  62.5  fim.  To  determine  whether  image 
vortices  introduced  by  the  symmetry  boundary  con¬ 
ditions  significantly  impact  these  results,  we  also  per¬ 
formed  calculations  on  larger  domains.  Results  on 
the  smaller  domain  depicted  in  Fig.  1  are  insensitive 
to  the  domain  size  for  the  first  10  ms  of  simulation 
time;  for  longer  times  there  are  noticeable  edge  ef¬ 
fects.  Where  we  report  longer  time  results,  we  use 
data  computed  on  domains  sufficiently  large  that  the 
edge  effects  are  not  manifest. 

The  inlet  fuel  for  the  simulations  consisted  of 
methane  uniformly  mixed  with  nitrogen-diluted  air 
(N2:02  =  6.28:1,  by  volume).  Five  stoichiometries 
were  considered:  4>  equal  to  0.76,  1.01,  1.10,  1.19, 
and  1.28.  In  Fig.  2,  we  present  snapshots  depicting 
the  mole  fraction  of  CH  for  each  stoichiometry.  At 
7.9  ms,  the  vortex  center  has  traveled  80%  of  the 
distance  initially  separating  the  flame  and  vortex  and 
has  interacted  with  the  flame  significantly.  The  ex¬ 
tent  of  the  simulated  decline  in  the  CIT  signal  is 
dependent  on  the  overall  flame  stoichiometry.  There 
is  no  noticeable  change  in  the  peak  value  of  CH 
along  the  flame  surface  for  </>  =  0.76  and  only  a  very 
small  decline  for  the  4>  =  1-01  case.  However,  at  (f> 
=  1.10,  the  decline  becomes  pronounced,  and  it 
steepens  with  further  increases  in  the  equivalence 
ratio. 

The  case  <j>  =  1.19  corresponds  to  the  stoichi¬ 
ometry  (and  nitrogen  dilution)  used  in  Nguyen  and 
Paul’s  experiment.  They  observe  a  two-  to  threefold 
decrease  in  the  averaged  CH  signal  before  it 
abruptly  disappears.  The  loss  of  signal  occurs 
throughout  the  flame  sheet  wrapped  around  the  vor¬ 
tex  and  occurs  about  7  ms  after  the  vortices  first 
begin  to  distort  the  flame.  Examination  of  our  nu¬ 
merical  results  indicates  that  CH  concentration  for 


this  case  exhibits  a  similar  rapid  decline,  leading 
eventually  to  a  nearly  10-fold  decrease.  This  simu¬ 
lated  rate  of  decrease  is  consistent  with  the  planar 
laser-induced  fluorescence  signal  descending  to  its 
noise  level  in  the  time  interval  observed.  Moreover, 
Fig.  2  shows  that  the  simulated  CH  concentration 
declines  precisely  in  the  region  identified  by  Nguyen 
and  Paul. 

A  more  detailed  examination  of  the  computational 
results  indicates  that  there  is  actually  a  qualitative 
shift  in  the  behavior  of  CH  between  <fi  =  1.10  and 
4>  =  1.19.  In  Fig.  3  we  show  time  histories  of  the 
peak  molar  concentration  of  each  species  along  the 
right  boundary  of  the  computational  domain  for  the 
two  cases.  In  the  rich  cases  (</>  £  1.19),  the  peak  in 
CH  declines  more  rapidly  than  the  peak  heat  re¬ 
lease;  the  opposite  is  true  in  the  lean  cases  (</>  £ 
1.10). 

In  this  region,  strain  effects  induced  by  the  flow 
field  components  tangential  to  the  flame  surface  as¬ 
sociated  with  the  impinging  vortex  result  in  a  strong, 
positive  flame  stretch.  Flame  curvature  caused  by 
flow  components  normal  to  the  flame  surface  result 
in  a  much  weaker  negative  stretch.  The  strain  and 
net  stretch  are  at  a  maximum  along  the  centerline 
(the  right  boundary  of  the  computational  domain) 
and  locally  decrease  as  we  move  along  the  flame  sur¬ 
face.  Classical  planar  strained  flame  analysis  (cf.  Ref. 
[36],  and  for  complex  chemistry,  see  Kee  et  al.  [37]) 
predicts  a  reduction  in  temperature  and  a  steepen¬ 
ing  of  gradients  with  increasing  strain  but  does  not 
predict  a  change  in  flame  structure.  The  experiments 
and  our  computations  corroborate  the  first  two  ef¬ 
fects  but  also  show  that  the  interaction  clearly  mod¬ 
ifies  the  flame  structure. 

We  note  that,  unlike  the  stretched  steady  flat- 
flame  case,  the  strong  circulation  pattern  around  the 
vortex  leads  to  strong  and  variable  tangential  flow 
along  the  flame  front.  As  a  part  of  this  process  the 
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Fig.  3.  Time  history  of  the  peak  values  of  molar  concentration,  normalized  to  their  respective  values  at  t  =  2  ms. 
Normalization  is  made  at  this  later  time  to  avoid  initial  transients  as  the  PREMIX  solution  adjusts  to  the  AMR  grid.  Note 
that  these  data  are  taken  from  simulations  on  very  large  1.2  X  5.6  cm  domains  to  avoid  edge  effects  at  later  times  in  the 
flame-vortex  interaction. 


circulation  of  the  vortex  entrains  material  from  the 
preheat  region  in  front  of  the  flame,  replacing  the 
equilibrium  distribution  of  combustion  intermedi¬ 
ates  upstream  of  the  primary  heat-release  zone  with 
the  inlet  mixture.  We  refer  to  this  specific  mechan¬ 
ical  process  as  one  that  “scours”  combustion  inter¬ 
mediates  from  the  front  of  the  flame.  The  net  effect 
of  this  process  depends  on  the  composition  of  the 
fluid  in  front  of  the  flame,  which  in  turn  depends  on 
the  equivalence  ratio. 

To  isolate  the  role  that  this  vortex-induced  scour¬ 
ing  has  on  flame  chemistry,  we  have  modeled  the 
process  in  one  dimension  with  source  terms  that  re¬ 
place  material  in  the  preheat  region  in  front  of  the 
flame  with  the  inlet  mixture.  This  one-dimensional 
model  thus  represents  the  chemical  effects  of  scour¬ 
ing  without  the  fluid  mechanical  effects  of  strain. 
Numerical  tests  with  this  one-dimensional  model  in¬ 
deed  show  that  for  the  richer  flame,  </>  =  1.19,  the 
scouring  effect  disrupts  the  chemical  pathways,  lead¬ 
ing  to  the  production  of  CH.  Conversely,  we  ob¬ 
served  essentially  no  change  for  the  lean  (f>  =  1.10 
case.  These  results  match  the  behavior  of  our  two- 
dimensional  study. 

Although  the  notion  of  vortex  scouring  in  combi¬ 
nation  with  effects  due  to  stretch  and  strain  describe 
the  mechanisms  by  which  the  vortex  modulates  the 
flame,  the  net  result  of  this  modification  can  be  quite 
complex.  It  is  compelling  to  conjecture  that  the 


dominant  effect  in  the  rich  case  is  the  removal  of  H2 
from  the  front  of  the  flame,  since  H2  profiles  in  the 
unperturbed  flame  extend  well  into  the  region  most 
affected  by  the  impinging  vortex.  However,  tests  of 
this  hypothesis  with  the  one-dimensional  scouring 
model  described  above  show  that  this  is  not  the  case. 
In  general,  the  effects  of  vortex  scouring  depend  on 
preferential  diffusion  effects  in  competition  with  re¬ 
actions  and  on  the  specific  reaction  pathways  that 
are  dominant  in  a  given  region  of  the  flame. 

The  flame  structures  were  also  analyzed  using  re¬ 
action  path  analysis.  The  reaction  pathways  are  those 
that  one  would  expect — namely,  the  Cj  pathway  to 
formyl  radical  and  an  enhancement  of  the  C2  path¬ 
ways  with  increasing  equivalence  ratio.  The  sensitiv¬ 
ity  analysis  showed  the  dominance  of  reactions  that 
are  most  important  in  controlling  the  size  of  the  rad¬ 
ical  pool,  and  in  many  cases  the  sensitivities  showed 
self-similarity,  indicating  that  only  a  few  variables 
control  the  state  of  the  system.  In  agreement  with 
Najm  et  ah,  we  find  that  the  formyl  radical,  HCO, 
is  a  reliable  marker  for  heat  release.  The  numerical 
results  presented  in  Fig.  3  and  similar  analyses  for 
other  stoichiometries  confirm  this  for  flames  inter¬ 
acting  with  vortices.  The  dramatic  decrease  of  CH 
in  rich  flames  is  linked  to  the  decrease  of  H  atom 
concentrations  and  triplet  methylene  and  decreases 
in  important  production  rates.  Interestingly,  some 
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Fig.  4.  CH30  profiles  along  the  vortex  centerline  for  </> 
=  1.19.  During  the  flame-vortex  interaction,  there  is  a 
marked  enhancement  of  CH30  concentration  in  low-tem¬ 
perature  regions. 

species  increase  in  concentration  in  low-tempera¬ 
ture  zones.  For  example,  for  (f>  =  1.19,  the  species 
CHsO,  C2H2,  H02,  CH3,  HCCOH,  CH2CO,  and 
C02  increase.  The  most  dramatic  example  occurs  for 
CH30,  as  shown  in  Fig.  4.  We  note  a  shift  of  the 
peak  concentration  to  lower  temperatures  as  the  vor¬ 
tex  couples  more  strongly  with  the  flame.  Sensitivity 
analysis  indicates  a  linkage  between  the  increase  in 
CH30  and  the  decline  in  CFI.  Another  noteworthy 
feature  is  the  increase  in  the  C2H2  concentration, 
which  occurs  as  a  result  of  the  flame-vortex  inter¬ 
actions.  This  prediction  is  potentially  significant  for 
soot  production. 

Conclusions 

Numerical  simulations  of  vortex-flame  interac¬ 
tions  have  revealed  a  strong  dependence  of  CH  be¬ 
havior  on  equivalence  ratio.  Analysis  of  this  behavior 
has  shown  that  the  vortex  not  only  stretches  and 
strains  the  flame  but  also  scours  material  from  the 
preheat  region  in  front  of  the  flame  and  at  early 
times  replaces  it  with  inlet  gases.  This  scouring  effect 
produces  the  dramatic  decline  in  CH  observed  ex¬ 
perimentally  by  Nguyen  and  Paul.  The  model  pre¬ 
dicts  a  measurable  increase  in  C2H2  as  a  result  of 
the  vortex  and  a  marked  increase  in  the  low-tem¬ 
perature  chemistry  activity. 

The  model  fails  to  reproduce  the  enhancement  of 
OH  observed  by  Nguyen  and  Paul.  One  possibility 


is  that,  with  our  current  limited  knowledge  of  the 
low-temperature  chemistry,  some  unknown  chem¬ 
istry  or  imprecise  knowledge  of  thermochemistry 
might  lead  to  the  experimentally  observed  rise  in 
OH  signal.  The  approximation  of  the  V-flame  con¬ 
figuration  using  a  flat  flame  and  enhanced  mixing  not 
represented  in  the  two-dimensional  flow  model  are 
other  possible  sources  of  discrepancy. 
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